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We consider the entropy associated to the phonons generated via the Hawking mech- 
anism in a sonic hole in a Bose-Einsten condensate. In a previous paper, we looked at 
the (l-l-l)-dimensional case both in the hydrodynamic limit and in the case when high- 
frequency dispersion is taken in account. Here, we extend the analysis, based on the 't 
Hooft brick wall model, by including transverse excitations. We show that they can cure 
the infrared divergence that appears in the (l-l-l)-dimensional case, by acting as an effec- 
tive mass for the phonons. In the hydrodynamic limit, where high-frequency dispersion 
is neglected, the ultraviolet divergence remains. On the contrary, in the dispersive case 
the entropy not only is finite, but it is completely fixed by the geometric parameters of 
the system. 
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1. Introduction 

Direct observations of the Hawking radiatioril' are extremely difficult, mainly be- 
cause h and are very small and this leads to a temperature of several orders of 
magnitude lower than the cosmic microwave background for a black hole of a solar 
mass. However, in condensed matter physics there are systems able to reproduce 
curved spacetime configurations, such that the speed of light is effectively replaced 
by the speed of sound waves. These are described by an equation of motion for- 
mally identical to the relativistic Klein-Gordon equation on a curved background-'^' 
In this context, an irrotational fluid accelerated from subsonic to supersonic speed 
generates a thermal flux of phonons with a spectrum analogous to the one asso- 
ciated to the Hawking radiation emitted by a gravitational black hole.'^'S These 
systems, known as acoustic black holes (or dumb/sonic holes), have attracted a 
great attention in the past years, especially because they might be soon realized in 
a laboratory. 

The analogy between astrophysical and acoustic black holes has been strength- 
ened thanks to the careful analysis of the correlation functions of the modes gen- 
erated via the Hawking mechanism near the acoustic horizon, both analyticallj'^'^ 
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and numericallyP The non-trivial structure of these correlation functions leads to 
ask what are the features of the entanglement entropy associated to these modes. 
For astrophysical black holes, it is well known that the area of the horizon is pro- 
portional to the thermodynamical entropyl^^J However, for acoustic black holes such 
a formula does not seem to exists. Nevertheless, it is possible to calculate the entan- 
glement entropy between the pairs of phonons created via the Hawking mechanism. 
In Ref. [HI we looked at this problem by considering the simplest system, namely 
a (l-l-l)-dimensional acoustic black hole in a Bose-Einstein Condensate (BEC). In 
this work, we used the so-called brick wall method, elaborated by 't Hooft several 
years ago and applied to gravitational black holes.'^^'^ At first sight, such a system 
seems trivial. In the realm of gravity, it is known that (l-l-l)-dimensional black holes 
are conformally invariant, and their entanglement entropy must be independent of 
the only characteristic scale of the black hole, namely its mass. In fact, in these 
systems the entropy is simply given by S" '-^ ln(L/e), where L and e are the infrared 
(IR) and the ultraviolet (UV) cutoff respectively.'^^'Sl Yot black holes in four dimen- 
sions, conformal invariance is broken, and the entanglement entropy turns out to 
scale as the area of the horizon, just like the thermal entropy. In the acoustic (1+1)- 
dimensional acoustic black holes studied in Ref. [TTJ we broke conformal invariance 
by introducing high frequency dispersion and we found non-trivial corrections to 
the entropy. 

Technically speaking, there is no rigorous proof that the entanglement entropy 
is the same as the brick wall entropy. However, there are strong evidences that this 
is the case, at least in the condensed matter systems considered here. To begin 
with, in Ref. [11] we showed that the two kinds of entropy are exactly the same in 
the conformal case, i.e. when no dispersion is taken in account. The introduction of 
high frequency dispersion leads to a correction to the leading term that is the same 
as the one found in spin-1/2 Heisenberg XX chains, using the "replica trick" , see e.g. 
Ref. [2TJ Finally, the results in Ref. [11] are in agreement with the ones found in Ref. 
1221 where the entanglement entropy was calculated using the density matrix p and 
the usual formula S = — Trplnp in a (l-l-l)-dimensional degenerate ideal Fermi gas 
with a sonic event horizon. Concerning astrophysical black holes, it has been shown 
that the two entropies are both proportional to the horizon area but they show a UV 
divergence that can be cured with the introduction of a cutoff around the Planck 
scale, see e.g. Ref.[T8l or a modification of the dispersion relation. Therefore it is 
possible that a proper regularization scheme makes entanglement entropy the same 
as the brick wall one. Such a suggestion is explored in Ref. [231 where the brick wall 
entropy is calculated for a Reissner-Nordstrom black hole and the divergences are 
regularized with the Pauli-Villar method. The authors found logarithmic corrections 
that are the same as the one-loop corrections to the Bekenstein entropy calculated 
with other methods. Later, is was realized that these corrections are the same as 
the one found with the "replica trick" whenever the background has vanishing Ricci 
scalar.E^ 
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In (l+l)-diniensional sonic holes in BEC, the UV cutoff is fixed by a specific 
dispersion function at high frequency that introduces a differential operator of order 
four in the mode equation for the phonon field outside the horizon, in analogy 
with certain gravitational models endowed with modified dispersion relations!^ 
In addition, when transverse phonon excitations are taken in account, they act 
as an effective mass that can remove IR divergences. Therefore, we expect that 
the entropy is completely determined by the microscopic parameters of these sonic 
holes, in opposition to the conformal case, where it depends only on the linear size 
and on an arbitrary cut-off. 

In the present work, we aim to prove this by extending the results of Ref. Illl to 
the case when transverse modes are excited. Typical laboratory BEC systems are 
made by a flow of ultracold atoms (e.g. Rubidium) modulated along a waveguide. If 
no transverse excitations are allowed, the system is essentially (l+l)-dimensional. 
However, transverse modes can be taken in account experimentally.'^^ 

The plan of the paper is the following. In Sec. 2 we lay down the basic equations 
that describe the acoustic black hole with the inclusion of transverse modes. In Sec. 
3 we calculate the entropy in the hydrodynamic limit in the presence of transverse 
excitations. In Sec. 4 we compute the entropy by including high frequency dispersion 
and we discuss our results in Sec. 5. 



2. The sonic black hole 

In the dilute gas approximation,'^ the BEC can be described by an operator ^ that 
obeys the equation 

ihdti' - f-|^v2 + v,y,t + 5*^*) * , (1) 

\ 2m J 

where m is the mass of the atoms, g is the non-linear atom-atom interaction con- 
stant, and Vext is the external trapping potential. The wave operator satisfies the 
canonical commutation relations [^(t, x), ^(i, x')] — S^{x~'x'). To study linear fluc- 
tuations, one substitutes 5* with ^o(l + 0) so that solves the Gross-Pitaevski 
equation 

ihdt'i'o = f-l-V" + Voxt + gn] ^-0 , (2) 



and the fluctuation (/) is governed by the Bogolubov-de Gennes equation 

thdt$ = -^ (v2 + 2^V ) $ + mc\$ + ^^), (3) 
2m \ Wo / 



where c = ^fgnjrn is the modulus of the speed of sound, and n = j^'op is the number 
density. We now focus on a one-dimensional flow of BEC that, conventionally, moves 
along the x-direction from right to left with constant velocity v — (— w,0, 0) and 
fixed n, while c— (— c(a;),0,0) smoothly varies in such a way that at the surface 
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X = we have c{x) = v. Conventionally, we choose c{x) < v{> v) for x < 0(> 0), 
thus we have a sonic horizon at cc = 0, with a supersonic region on the left and a 
subsonic one on the right. With this choice, v and c(x) are positive numbers. 

It is possible to adjust c{x) provided one allows for a varying g, such that the 
quantity gn + Vcxt remains constant in space.^^In this way, Eq. ([2]) admits the plane- 
wave solution ^0 ~ \/nexjp(iko ■ x — iujQt), where /cq is related to the condensate 
velocity through the relation v — hko / m. To study the dynamics of i/i, it is convenient 
to consider the so-called density-phase representation defined by 




along the lines of Ref. 5 By choosing a one-dimensional flowing condensate, we 
reduce the acoustic metric to the fornP 

ds'^ = — ^ l-ic^ix) - v^)dt^ - 2vdxdt + dx^ + dy'^ + dz^] , (5) 
mc[x) 

which represents an acoustic black hole in (1-fl) dimensions, implemented by two 
asymptotically flat directions. In fact, this metric has an event horizon located where 
c{x) = V and its structure is the same of the Painleve-Gullstrand line element,'^^ 
up to the conformal factor n/(mc). Normally, one discards transverse directions by 
considering only a;-dependent phonon fields and by neglecting the (y, z) metric 
components. In fact, the physical model is an elongated cloud of ultracold atoms 
trapped by means of magnetic fields and moving with constant speed along the 
a;-axis and transverse vibrations are usually ignored.!^ Here we propose to take in 
account these additional transverse degrees of freedom by assuming that the phonon 
field has the expression 

9 = 0{x) ex:p{ikyy + ikzz) , (6) 
h — h(x) exp{ikyy + ik^z) . (7) 

If we now expand the operators and fi into sums of creation and annihilation 
operators, i. e. in the form 

n = ^(an-l-h.c.), ^ ^(a ^ + h.c.) , (8) 

we then find the two coupled equations 

o^{dt + vd,)h + -die - -^e , (9) 

where A"^ = fc^ + fc^. The function ^ = h/{mc) is known as the healing length, which 
sets roughly the scale at which the model breaks down.^^Sl In these expressions, we 
see that the effects of the transverse modes appear as an effective mass, that can 
be quantized in the case when the (y, z) range is bounded or periodic. Note that 
we would have obtained these very same equation by applying the Klein-Gordon 
operator to a massive scalar field propagating only on the (t, x) sector of the acoustic 
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metric as in Ref. 12, Therefore, the term A^^ can be correctly considered as an 
effective mass for the dimensionaUy reduced field and the system is equivalent to a 
massive scalar fields propagating on a (l+l)-dimensional black hole. In particular, 
the wave equation has no angular quantum numbers, as it happens for (3+1)- 
dimensional black holes and in the calculation of the entropy, we expect that the 
contribution of the extra-dimensions is completely encoded in the mass term. This 
is also supported by the fact that the horizon has no area, it is just a point located 
at a; = 0. 

In what follows, we will consider phononic modes with wavelength much smaller 
than ^. Note that, in our settings, ^ is not constant, as it depends on the (varying) 
speed of sound. The limit ^ = corresponds to the so-called hydrodynamic limit. 
Within this regime, one can show that the two equations above can be combined 
together into 



n{x) 



{dt + vdx) 



1 



c(x)- 



e{x) 



{dt+vd^)e{x) 



(10) 
(11) 



3. Brick wall and entanglement entropy in the hydrodynamic limit 

In this section, we study the equation ([T0|) and use the 't Hooft brick-wall method 
to compute the entropy.GSl To begin with, we assume that the field 9 is stationary, 
and wc set 6 = 9{x, A) exp{iujt). Then, Eq. ([TUl) can be written in the form 

e"{x, A) + 2A{x)9\x, A) + B{x, X)e{x, A) = , (12) 

where 



A{x) 



iojv + v 



1 c' 



B{x,X) = 



2iujv- 



(13) 



— — 
For notational simplicity, from now on the x-dependence of c(x) will be understood. 
To apply the WKB formula, we make the substitution 

1 



e{x,\) 



: exp 



^ J f{x,X)dx 



where 

f{x,X)^gix,X) + ihA{x). 
In this way, Eq. (|12|) assumes the WKB form 



where 



3g;^ 
452 

Pix) = 
Vix,X) - 



El 



In" , , V(x, X) 
-^ + P(.) + ^^0 



(14) 



(15) 



(16) 



w2(3c2-2t;2) fc 



(c^ — i;2)2 



(c^ — v"^) c 



2x2 



CJ^A' 



(c^ — W^)^ (c^ — t|2) 



(17) 
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The mode density number is obtained by integrating f(x, A) over the near-horizon 
region. If we assume that the transverse momenta ky and kz are quantized, we also 
need to sum over aU their possible values, in analogy with the sum over the angular 
momenta in the gravitational (3 + l)-dimensional case. 

The WKB analysis proceeds by substituting g = gQ + h?g2 + ti^g^ + . . . in the 
above equation and by collecting the terms with the same power of h. The lowest 
order is 

g^{x,\)^±^V{x,\) , (18) 

so the lowest WKB order solution to Eq. reads 

/o(x,A) = ±^/V{x,\)+ihA{x) . (19) 

We find convenient to sum first over the transverse momenta. By turning the 
infinite sum into an integral we have 



foix) = r!'^2 (2^')^ I dkydk.^l-Z^{kl + kl)+zhA{x), (20) 



c2 - (27r)2 

where we set 



c 



2 



^ = — (21) 



and a / (27r)2 is the surface spanned by the transverse oscillations. The sign has been 
fixed so that we recover the known case for ky,kz =0. By integrating in the region 
where the square root is real, we find 

= 2^^ + 2W^- 

To obtain this formula, we used the approximation c = v+kx, where k = {dc/ dx)x=o 
is the analog of the surface gravity of a black hole. This is justified by the hypothesis 
that we want to integrate the mode number in the near-horizon region, where c{x) 
is linearized. The same approximation is usually applied to the (3 + l)-dimensional 
gravitational black hole. 

We now compute the entanglement entropy by integrating over the number of 
modes with energy E = fno living in a segment with one endpoint placed at an 
arbitrarily small distance x — e from the horizon and the other at x = L, where L 
can be identified with the length of the region where the near-horizon approximation 
is valid. It is easy to show that i = u/k is a good estimate for the length of this 
region, provided k is a smooth function of x. We find 
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We now recall that the free energy and the entropy associated to massless spin-O 
particles are given respectively by 

p n{E)dE 

where /3 = (/cbT)^^ = 2tt /{Hk) is the inverse of the (black hole) temperature. With 
these formulae, we find that the entropy reads 

We see that the transverse modes contributes to the entropy with a geometric term 
that still diverges in the UV (e -> 0). If we allow for a transverse section with a 
radius of the order of i, a — vrL^, we see that this contribution is negligible for 
L ^ e. Had we not taken in account these, we would have founcP^ 

Sx^^^Un(^^y (27) 

One might be worried because the expression (|26l) does not converge to this result 
for cr — 0. This is not so surprising as the limit A — )■ oo and the integral in 
Eq. ([20l) do not commute. The same situation occurs in the usual 't Hooft brick 
wall model: if one does not integrate over the angular momentum by choosing, for 
example, i = the result is radically different and one does not even recover the 
proportionality between horizon area and entropy. The conclusion of this section is 
that the transverse mode contributes to the entropy but are not able to cure the 
ultraviolet divergence in the hydrodynamic limit. The situation is quite different 
when dispersion is taken in account, as we show in the next section. 



4. Dispersive case 

We now consider the dispersive case, and evaluate the contribution to the entropy 
given by the transverse excitations. In order to do so, it is more convenient to write 
the Eq. ([3]) by using the expansion (p = ^^ [a^ + a^t^*]. By assuming, as above, that 
4> = 4>(x) exp{iijjt) e'xjpliky'y + ik^z) and that (p ~ '~p(x) exp(icji) e:sq){ikyy + ik^z), we 
obtain the two coupled equations 



i-cu + zvd,)^ + ^dl4> - ^0 - -(^ + ^) = 0, 



{uj - ivd.,)^ + - - |('^ + ^) - 0. (28) 

By following Ref. Ill[ we decouple these equations in momentum space. By defining 

r ,irj.p-ipx 

V±(x) = ^{x) ± ^{x), i,^{p) - / _^V^(a:), (29) 
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we find 



(P) 



2mX^{uj — pv) ~ 



where 



V{p, A) - 



h{p^x^ + 1] 
\'^{uj-pvf _ 



p^\^ + 1 

A2 



(30) 
(31) 

(32) 



Eq. ([50)1 correctly reduces to the one found in Ref. [TT] when A ^ oo and we can 
solve it in the near-horizon approximation, i.e. when c = v + indp. In analogy to 
what we did in the hydrodynamic limit, we use a WKB procedure and replace the 
expression 

1 



Vf¥) 



exp 



' dp ( f{p) + ^ 



in Eq. ^ to find 



3 /ry r 



(33) 



(34) 



where the prime here denotes differentiation with respect to p. Note that the term 
vh/ K is not present in the momentum representation used in Ref. [11] This is due to 
the fact that here we use the (</>, (p) representation while in the hydrodynamic limit 
of Ref. [TT] we use the density-phase representation (ni, 6*1). 
By keeping the lowest order only, we find the solution 



1 ± 



A2p2 



AV 



1 



vp J 4 



AV + 1 
A2 



(35) 



where ^0 = ^(^^ ~ 0) is the healing length at the horizon. By taking the hydro- 
dynamic limit ^0 0, and by computing the number of modes, we see that, in 
order to recover Eq. (|23p we must consider only the sign a choice that will be 
understood in the following. For ^0 7^ Oj we can check that, in the limit A 00, 
this equation becomes the same as Eq. (21) of Ref. [TT]as expected. However, the 
role of A turns out to be crucial. In fact, we note that Eq. ([35]) is no longer singular 
in the limit p — )• 0. This means that in the computation of the number of modes we 
do not have to worry about the low-momentum modes, which had to be eliminated 
by a IR cut-off in the calculations of Ref. [TT] Therefore, we see that the effect of 
the transverse modes is to naturally regularize the infrared divergence. 

The number of modes with a given energy E — huj follows from the formula 

1 /-P.. 



i{E) = 



nh 



A 



f{p, \)dp, 



(36) 



where the interval of integration is determined by the condition that the square root 
in /(p. A) is real. Before integrating over p, we need to calculate the sum over the 
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possible values of A. As we did in the hydrodynamic limit, we turn the sum over A 
into an integral so that 



AE) = — / 



dp 



1 



Sir 



where 



H = 



(- 

V V 



P 



,(37) 



(38) 



and Ep is the incomplete elliptic integral of first kind. The range of integration 
[0,p] coincide with H G [l,c»[. Exact integration is not possible, but, by noting 
that the integral is dominated by large values of H, approximations can be made 
and the mode density number can be written as 



/ r^N E La \ - 



hv 



(39) 



where, again, we set L = v/k as the size of the near-horizon region, and aj are 
numerical coefficients. The index j run over integers greater than 4 except for the 
first two terms of the sum, for which j ~ 5/2 and j = 3. In any case the terms of 
the sum are subdominant in the calculation of the entropy. In fact, by using the 
formulae and we find that the entropy has the form 



S = 



1 a 



^5/2 



5/2 



as 



3=5 



(40) 



where the a^-'s differ from the a^-'s by unimportant numerical factors. The important 
result instead is that the leading term is a constant and the logarithmic dependence 
in Eq. (l26t has disappeared. The correction terms are all small because of the 
condition that the healing length at the horizon is much smaller than the size of 
the near-horizon region L. Again, if we allow for a transverse section with radius of 
the order of the near horizon region, we see that all the term in the sum have the 
form (Co/i)'" with A; > 0. 



5. Discussion 

Let us briefly summarize the main results. In the hydrodynamic limit, we computed 
the entanglement entropy by taking in account transverse excitations. As expected, 
these act as an infrared cutoff, and the corrections to the usual logarithmic term 
(l/6)ln(L/e) have the form a/Le. The presence of the transverse modes does not 
help in curing the UV divergence encoded in the limit e — > and that represents 
the effect of modes boundlessly piling up near the horizon, as a result of the lack of 
any ultraviolet completion of the field theory. 

In this perspective, we find very interesting the results of Sec. 4, where dispersion 
effects are taken in account. For a BEC the dispersion relation is modified at high 
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frequency according to the functiorP 



{uj — vpY 




) 



(41) 



In Ref. [TT]we have found that the modified dispersion relation acts as an ultraviolet 
cutoff and the entanglement entropy receives a correction that depends on the heal- 
ing length. The leading term is however still arbitrary, as it depends on the largest 
wavelength allowed in the system. This reflects the lack of a physical infrared cutoff 
that can be provided for by taking in account the effects of transverse modes. As we 
have seen in Sec. 4, the density of modes turns out to be finite in the near-horizon 
region, thus revealing a surprising IR/UV mixing: the infrared cutoff A prevents the 
divergence of the density of modes near the horizon. This has a dramatic conse- 
quence as the leading term of the entropy becomes a constant and the corrections 
are uniquely fixed by geometric parameters of the system: the healing length ^q, 
the extension of the near-horizon region L, and the cross section <j. One challenging 
issue is to verify that the entanglement entropy found here is the same as the one 
calculated by other methods. We hope to report soon on this topic. 
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